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Electron Spin Resonance (ESR) describes a suite of techniques for characterising electronic sys¬ 
tems, with applications in physics, materials science, chemistry, and biology. However, the require¬ 
ment for large electron spin ensembles in conventional ESR techniques limits their spatial resolution. 

Here we present a method for measuring the ESR spectrum of nanoscale electronic environments by 
measuring the relaxation time (Ti) of an optically addressed single-spin probe as it is systematically 
tuned into resonance with the target electronic system. As a proof of concept we extract the spec¬ 
tral distribution for the PI electronic spin bath in diamond using an ensemble of nitrogen-vacancy 
centres, and demonstrate excellent agreement with theoretical expectations. As the response of 
each NV spin in this experiment is dominated by a single PI spin at a mean distance of 2.7 nm, the 
extension of this all-optical technique to the single NV case will enable nanoscale ESR spectroscopy 
of atomic and molecular spin systems. 


I. INTRODUCTION 

Techniques to detect electron spin resonance (ESR) 
have long been used to study materials and systems con¬ 
taining unpaired electron spins, such as metal complexes 
and organic radicals. From an operational viewpoint, 
the low spin density and high decay rates of such radi¬ 
cals in-situ, together with the limited sensitivity of ESR 
detection systems, often makes the task of obtaining an 
ESR signal above the detection limit difficult. This is fur¬ 
ther complicated by the fact that a major component of 
any biological tissue is water, which has a strong electric- 
dipole mediated absorbtion band in the same region of 
the electromagnetic spectrum (microwave) as the signals 
emitted in ESR experiments^. As such, there is a great 
need for a highly sensitive, highly localised technique that 
may be used to obtain ESR spectra from unpaired elec¬ 
tron systems, without the need for microwave control of 
the sample. 

The nitrogen-vacancy (NV) centre in diamond is a 
promising candidate for such an ideal nanoscale ESR 
probe. Much attention has been focused on using 
measurements of NV spin energy sublevels^— or de¬ 
phasing rates (1/T 2 )£~— to characterise dynamic pro¬ 
cesses occurring in external magnetic environments. 
These protocols have been shown to have remarkable 
sensitivity to frequencies in the kHz to MHz range, 
and are thus well suited to characterising nuclear spin 
environments^—. However, to achieve the desired sen¬ 
sitivity to the more rapidly fluctuating fields associated 
with electron spin environments, and more importantly, 
the ability to be frequency-selective, complex pulse se¬ 
quences are necessary^—. 

Here we focus on on a new, and potentially simpler, 
technique for extracting the ESR spectrum based on the 
longitudinal decoherence (or relaxation) time, Ti, of the 
NV centre. Given the GHz transition frequencies of the 


NV spin in its orbital ground state, owing to its zero- 
field spitting of D = 2.87 GHz (Fig. 1), paramagnetic 
spins can enhance the longitudinal decay rate of the NV 
spin, enabling one to detect their presence by monitor¬ 
ing changes in T)— — . Magnetic sensing using relaxation 
times can be much more sensitive for detection because 
Ti times of an NV centre can be up to three orders of 
magnitude longer than spin-echo based T 2 timesiS. 

In the presence of a controlled external magnetic field, 
B 0 , the transition energies of the NV probe, wnv = 
27 tD ± 7H0 (where 7 = 1.76 x 10 11 rads _1 T _1 is the 
electron gyromagnetic ratio), can be brought into reso¬ 
nance with environmental spin energies, we, via the Zee- 
man effect (Fig. [I]), and the corresponding change in Ti 
may be measured. We establish a general and robust 
method to determine the environmental spectral distri¬ 
bution 5 (we) by measuring the relaxation time of the 
NV probe as a function of the external field strength, 
Ti(Bq). Our method is tested experimentally by measur¬ 
ing the ESR spectrum of the substitutional nitrogen (PI 
centres) donor electron spins in type-lb diamond over 
the range 0-200 MHz. Comparison with theoretical ex¬ 
pectations for this known system confirms the validity of 
this general approach for measuring the ESR spectrum 
using the NV centre spin-probe in a range of nanoscale 
applications. 

The nitrogen-vacancy centre (see Refs. [2^ and (26| for 
extensive reviews) is a point defect in a diamond lat¬ 
tice comprised of a substitutional atomic nitrogen im¬ 
purity and an adjacent vacancy (Fig.|T](a)). The energy 
level scheme of the C3„-symmetric NV system (Fig.Q](b)) 
consists of ground ( 3 A 2 ), excited ( 3 E) and singlet elec¬ 
tronic states. The ground state spin-1 manifold has three 
spin sub-levels (| 0), |±l)), which at zero field are split 
by 2.87 GHz. Since the magnetic sublevels have different 
fluorescence intensities when illuminated with laser light, 
we can achieve spin-state readout opticall y^ 23 . The 
degeneracy between the |±l) states may be lifted with 
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FIG. 1: (a) The nitrogen-vacancy (NV) centre point defect 
in a diamond lattice, comprised of a substitutional atomic ni¬ 
trogen impurity (N) and an adjacent crystallographic vacancy 
(V). (b) The NV ground state spin sublevels are separated by 
D = 2.87 GHz. Upon optical excitation at 532 nm, the pop¬ 
ulation of the jo) state may be readout by monitoring the 
intensity of the emitted red light, (c) Transition frequencies 
of the NV centre, and that of a single electron as a function 
of external magnetic field strength, Bo. At Bo = 512 G the 
|0) <->■ | — 1) transition of the NV is resonant with the electron 
transition, allowing them to exchange energy, (d) Illustration 
of the way in which the NV relaxation filter functions, T^ 
and Tf (having width r 2 = I/T 2 ), probes specific regions of 
the environmental spectral density. One of the NV transi¬ 
tions is significantly detuned from the frequencies exhibited 
by the environment, and hence cannot be excited. The NV 
spin thus behaves similar to a spin-i qubit and the steady 
state of the NV spin will be an equal mixture of only the two 
states associated with the excited transition, with the third 
state remaining unpopulated. 


the application of an external field, with a corresponding 
separation of 2xl7.6xl0 6 rads~ 1 G _1 , or 5.60 MHz G -1 , 
permitting all three states to be accessible via microwave 
control, however the |±l) states are not directly distin¬ 
guishable from one another via optical means. By isolat¬ 
ing either the | 0) «A |+l) or | 0) 0 | — l) transitions, the 
NV spin system constitutes a controllable, addressable 
spin qubit. 

The large zero-field splitting is fortuitous in the present 
context in that transitions in the ground-state spin triplet 
manifold are far detuned from the ^MHz couplings to ni¬ 
trogen donor electron and 13 C nuclear paramagnetic im¬ 
purities within the crystal, leaving NV transitions unable 
to be excited unless brought into resonance using an axial 
magnetic field ( Bq ps 512 G, for an electron spin environ¬ 
ment). The weak spin-orbit coupling to crystal phonons 
(and low phonon occupancy, owing to the large Debye 
temperature of diamond) leads to longitudinal relaxation 
of the spin state on timescales of roughly Xj ~ 1 — 10 ms 
at room temperature. On the other hand, the trans¬ 
verse relaxation time (T 2 ) of the NV spin is the result 


of dipole-dipole coupling to other spin impurities in the 
diamond crystal, which occur on timescales of 0.1-1 fis 
in type-lb diamond. Sensing based on the much longer 
relaxation times can thus be significantly more sensitive 
than those based on T 2 . Shot-noise statistics dictate that 
the minimum magnetic field detectable over a total in¬ 
terrogation time T = Nt, comprised of N cycles of dark 
time r, is proportional to 5 m i n oc 1/y/Tr. As the dark 
time, r, is limited by Tj in relaxation based protocols, 
the latter offer an improvement in sensitivity by a factor 
over the former for a fixed interrogation time, 
T. Hence, the relatively long relaxation time of the NV 
ground state and inherent sensitivity to GHz frequencies, 
together with room temperature operation and optical 
readout, make it an ideal system for the ESR spectral 
mapping protocol discussed in this work. 


II. RESULTS 

A. Relaxation of an NV spin in an arbitrary 
magnetic environment 

At zero field, both NV transitions, | 0) -A |+l) and 
I °) ^ |-1>, occur at approximately 2.87 GHz. These 
transitions are then split symmetrically by the presence 
of an axial magnetic field and any axial NV-environment 
couplings, and move to higher and lower frequencies, re¬ 
spectively, with the resulting shift being directly propor¬ 
tional to the strength of the field (Fig.QJc)). When the 
transition frequencies of a spin-1 NV system are brought 
into resonance with those of a particular environmental 
transition, the two will be able to exchange energy via 
their magnetic dipole interaction. As the shifts induced 
by the environment, typically of order MHz, are con¬ 
stantly fluctuating, they act to broaden the NV’s transi¬ 
tion frequencies, as characterised by their resulting inho¬ 
mogeneous linewidth, or transverse spin relaxation rate, 
T 2 (which may be either T 2 = 1 /T% or 1/T 2 depending 
on the whether an ensemble of NVs is used, or just a 
single NV spin). This results in a resonantly enhanced 
response of the NV spin relaxation rates (T 1 * 1 ) to envi¬ 
ronmental fields fluctuating at (or within T 2 of) the NV 
frequency (Fig.[Tt d)). 

To compute the response of the NV spin to such an en¬ 
vironment, we must account for two dominant processes: 
energy exchange between the NV and the environment 
(effective coupling rate B), which changes the popula¬ 
tion of the magnetic sublevels of the NV ground state; 
and the destruction of the phase coherence between these 
sublevels (depliasing, occurring at a rate of T 2 ). Depend¬ 
ing on the relative strengths of these processes (that is, 
how many energy exchanges may occur before the NV 
spin is dephased), the NV spin can exhibit diverse be¬ 
haviour (Fig. [2]). However, as there are typically more 
sources of dephasing than energy exchange (with the lat¬ 
ter effect further decreasing when the NV spin and envi¬ 
ronment are away from resonance), we have T 2 B for 
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Dimensionless time, t!T 2 [AU] 

FIG. 2: Bloch vector behaviour at constant dephasing 

rate (F 2 ), for cases of under-damped (B > T 2 / 2 \/ 2 , red), 
critically damped (B = F 2 / 2 \/ 2 , green), and over-damped 
(B < T 2 / 2 V 2 , blue) NV-Environment couplings, with the +2 
axis representing the |0) state of the NV spin, and the —2 
axis representing the | — 1) state, (a) Example realisations 
of Bloch sphere trajectories, (b) Average Bloch vector tra¬ 
jectories. (c) Population of the NV |0) state of the NV spin 
(P 0 (t) — ^ + 52 (f)). In the under-damped case the Bloch 
vector is able to rotate appreciably before its transverse pro¬ 
jection decays. In the over-damped case the transverse pro¬ 
jection of the Bloch vector is always pulled back to the z axis 
before it can appreciably rotate. Most practical cases reside 
in the latter regime, due to the intrinsic dephasing arising 
from strong coupling of the NV to paramagnetic defects in 
the diamond crystal. 


the cases considered in this work. 

Under this regime, if the NV is initially polarised in its 
| 0 ) state, and only the | 0 ) | — 1 ) transition is being ex¬ 

cited, the subsequent population at time t when coupled 
to a spin system of transition energy we is given by (see 
Supplementary Material for details) 

Po(uE,B 0 ,B,t) = i + iexp (-2B 2 ^ 2( ^ 2 +ri ),(l) 

where S = wnv — we is the difference in transition fre¬ 
quencies between the NV (wnv) and the environment 
(we)- 

Typical electron spin environments will exhibit a dis¬ 
tribution of coupling strengths ( B ) and frequencies (we), 
denoted Pb{B) and S (cue, -Bo) of roughly the MHz-GHz 
regime ^ 20 ' 29 , which must be taken into account to deter¬ 
mine the full response of an NV ensemble. This response 


is given by 

(Po{Bo,t)} = J G(we, B 0 ,t) S(cue, B 0 ) dcuE, 

where we identify the function 

G(cu E ,-B 0 ,f) = J P 0 (co E ,B 0 ,B,t)P B {B)dB, (2) 

as an environmental noise filter with a Lorentzian point- 
spread function centred on the NV transition frequency 
(Eqn.Q]). The NV-relaxation filter, G, is tunable via 
the Zeeman interaction, meaning that it may be directly 
tuned to specific parts of the environment’s spectral dis¬ 
tribution, S(u} E ,Bq), by choosing the strength of the ex¬ 
ternal magnetic field, Bq. 

If the environment is comprised of some distribution of 
single electrons, their transition frequencies, we (includ¬ 
ing spectral features such as hyperfine interactions, but 
centred about the 7 -B 0 Zeeman shift), may be brought 
into resonance with those of the | 0 ) -H- | — l) transition 
of the NV (i.e., cu_ ~ 2ttD — 7 S 0 ) by choosing a magnetic 
field strength Bo such that Bq = irD/'y ~ 512 G. In other 
cases, the environment may contain spin -1 or greater sys¬ 
tems possessing their own zero-field splitting (See Ref.i20l 
for the case of spin-| Gd spins coupled to individual NV 
centres). If the environmental zero-field splittings are 
greater than that of the NV, the | 0) <-> |+l) transition 
of the NV will need to be utilised to ensure that the re¬ 
spective energies are brought into resonance. 

In the following, we discuss how we may take advantage 
of the relaxation filter to reconstruct the spectral density, 
S(oj Ei Bo) of an arbitrary environment. 

B. Reconstruction of the environmental spectral 
density 

As noted above, the region of the spectral density sam¬ 
pled by the NV filter functions may be tuned by con¬ 
trolling the strength of the static external field. This 
suggests that, by sweeping the filter function across the 
entire spectrum, we should in principle be able to recon¬ 
struct it by measuring the relaxation rate of the NV spin 
for an appropriate range of external field strengths. 

We denote an arbitrary given noise spectrum at zero 
field by So (cue) = S’(cue, 0). The distribution at some 
finite external field B 0 is then So (cue — jB 0 ). For most 
cases of practical interest, we assume the shape of the 
distribution does not change with 7 .Bo (see Supplemen¬ 
tary Materials for more details), although this case can 
be handled by extension. Furthermore, we also assume 
that one of the NV transitions is sufficiently detuned that 
it is not sensitive to the environment, making the overlap 
with the spectrum insignificant. Even if this is not true, 
the detuned filter function will translate with 7 .Bo at the 
same rate as So (cue — 7 -B 0 ) (Fig.QJd)), and thus produce 
a constant shift in the overall measurement that does not 
change with Bq, which may be later subtracted. 
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FIG. 3: Sequence (a) By controlling the external field strength, Bo, the NV filter function, G may be tuned to filter specific 
regions of the spectral density, S. (b) The resulting measured signal, M(Q, o), is the convolution of the NV filter function, G, 
and the spectral density of the environment, S. (c) Given that G is known, the spectral density, S, may be reconstructed by 
deconvolving S and G from M. 


In principle, the filter function, G is fully defined by the 
hyperfine shift and free induction decay (FID) time (TJ) 
of the NV spin (see Supplementary Material for details), 
which may be characterised by performing a 7r/2 —t—7t/2 
FID (or Ramsey) measurement. The full G(we) function 
may also be mapped to arbitrary precision by measuring 
the response of the NV longitudinal relaxation rate to an 
applied oscillatory magnetic field of frequency we, pro¬ 
vided the microwave field strength, B (playing the same 
role as the effective coupling strength of the NV spin to a 
proximate spin system), is such that B <C 1/T 2 * (to avoid 
microwave broadening of G ). 

The measured response of the NV relaxation rate, 
to So for some external field strength, B 0 
(wo = 7 -Bo), is then given by (for brevity, we put 
M(ivo,t) = (Bo(t, B 0 )) meag ) 

/ OO 

S 0 (w E - wo)G(we - 27 tD + to 0 , t) dw E - 

-OO 

( 3 ) 


By introducing the frequency-space variable D = w E -wo, 
and the parameter, Do = 27 tD — 2wo, and making use of 
the symmetry properties of the Lorentzian function, we 
may write this integral as a Fourier-space convolution, 

/ OO 

So(D)G(D 0 — D, t ) dD 

-OO 

= (S 0 *G) t (D 0 ). (4) 

Given that the filter function is known, the spectral den¬ 
sity may thus be reconstructed using an appropriate de- 
convolution algorithm (see Supplementary Material for 
details). 

In the following section, we provide an experimental 
example of this protocol applied to the PI electron spin 
bath of type-lb diamond. 


C. Experimental demonstration 


In order to demonstrate the reconstruction of an en¬ 
vironmental spectrum using this technique, we measured 
the relaxation rate of NV centres subject to a bath of sub¬ 
stitutional nitrogen donor (or PI centre) electron spins. 
The sample chosen was a synthetic diamond grown with 
high-pressure high-temperature (HPHT) method, with 
~50ppm substitutional nitrogen and 1-10 ppm NV - con¬ 
centration. NV centers are oriented along all four symme¬ 
try axes, and a solenoid was used to apply an adjustable 
magnetic field along one of these orientations. 

After initializing NV spins in the |0) state via optical 
pumping (or the | —1) state with an additional microwave 
7T pulse) we measured the fluorescence following a vari¬ 
able amount of ‘dark time’ (for which the pump laser was 
off). NVs not aligned with the external field and other 
defects contribute to the diamond fluorescence, which 
may decay at rates different from that of the aligned 
NVs. In order to remove this effect, we employed a form 
of common-mode rejection (described in Ref.[30i). After 
initializing to either the |0) or | —1) state, we measured 
the fluorescence contrast after each dark time, and sub¬ 
tracted the results for each initial state. The resulting 
decay of the fluorescence difference indicates the loss of 
NV population from the initial state, and the magnetic 
noise these populations are sensitive to. We note here 
that for cases in which the application of microwaves may 
be impractical, this technique may be performed in an 
all optical manner, although this approach will also re¬ 
quire the characterisation of other sources of fluorescence 
within the sample. 

The PI centre, the target of our demonstration, is a 
substitutional defect of a nitrogen atom in place of a car¬ 
bon atom in a diamond lattice. The single unbonded 
electron able to reside along any one of the four crystal¬ 
lographic bond axes, giving rise to four possible hyper- 
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FIG. 4: (a) Theoretical plot showing the transition frequencies of the NV centre and PI electron spins vs the strength of an 
external field aligned along the (111) axis. Energy exchange between an NV spin and a nearby PI spin is achieved when 
two transition frequencies approach resonance. Broadening of these lines is caused by interactions with other PI spins and 
spin impurities within the diamond crystal, (b) Zoomed plot of that in (a), highlighting the liyperhne structure and the 
corresponding spin conserving and non-spin conserving transitions of the PI centre, c) & d) Plots of NV centre contrast vs 
dark time at 509.5 G and 481.3 G respectively. In the case of the former, NV and PI electron energies are on resonance, 
resulting in a comparatively rapid decay. In the case of the later, NV and PI transition energies are too far detuned to facilitate 
a dipole-mediated resonant energy exchange, meaning that the depolarisation of the NV spin ensemble is dominated by 
interactions with crystal phonons, e) Measurements of curves such as those shown in c) and d) for 500 magnetic field strengths 
between 480 G and 540 G. f) Plot of the external field-dependent relaxation rate, Pi(Bo), extracted from the data shown in 
e). Correspondences of the observed features with the PI spin transitions in b) are indicated explicitly, g) Application of the 
deconvolution procedure to the data in f) yields the spectral density of the spin bath environment surrounding the NV centres in 
the diamond sample (blue), showing good agreement with theoretical predictions (red - see Supplementary Material for details). 
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fine coupling configurations with the nitrogen nucleus. 
For field strengths above ~ 100 G, the quantisation axis 
of the PI electron spin is set by its Zeeman interaction, 
effectively reducing the number of possible PI species 
from four to two. Hence, there exists a 25% chance 
that the delocalisation axis of the PI centre is aligned 
with the NV axis, giving an axial hyperfine coupling of 
114 MHz, and a 75% chance that the delocalisation axis 
is arccos (— 5 ) ~ 109° to the NV axis, producing an ax¬ 
ial hyperfine coupling os 86 MHz (the spin properties of 
this defect are discussed in detail in the Supplementary 
Material). Figures. 0(a) & (b) show the overlap of the 
transition frequencies associated with both the NV spin 
and the PI centre. 

Measurements of the longitudinal spin relaxation of the 
NV ensemble, (Po(t, -Bo)) meas > were taken at 500 differ¬ 
ent external magnetic field strengths between 480 G and 
540 G (Figs.[H(c) to (e)). Before we deconvolve the data 
set to determine the spectral density, we extract and in¬ 
spect the magnetic field-dependent component of the NV 
relaxation rate due to spin-spin relaxation with the en¬ 
vironment. Accordingly, the data was fitted using the 
function given by 

P 0 fit (H 0 ,t) = exp(— y/TjBo)i-Rt) , (5) 

where R ~ 360 Hz is the component of the relax¬ 
ation due to environmental phonons (see Supplemen¬ 
tary Material for the derivation of this fitting form). 
The resulting spin-spin relaxation rates are plotted in 
Fig-IH(f)- As expected, the measurements show that 
the NV relaxation rate increases when D — 2ujq/2t: = 
0, ±86 MHZ, ±114 MHz, or B 0 = 490,495,510,525, 
and 530 G, corresponding to the conditions under which 
the NV spin may exchange magnetisation with PI 
electron spins. Other features are evident at Bq = 
499, 502, 517 and 520 G, which correspond to processes 
that change the magnetisation of the NV spin and the P1 
nuclear spin, but leave the PI electron spin unchanged (as 
seen previously in Ref. [31). These transitions are not en¬ 
ergy conserving and are thus comparatively weak, leaving 
them partially obscured by the large spin-phonon relax¬ 
ation effect of the NV spin, resulting from two-phonon 
Orbach^, and two-phonon^ Raman processes. 

To obtain the spectral distribution, we carry out the 
deconvolution over the entire set of data, (Po(t, Bo)) , 
and plot the spectral distribution in frequency space, 
S(f), in Fig.H](g), demonstrating good agreement with 
the corresponding theoretical expectations (see Supple¬ 
mentary Material for details). The deconvolution to the 
spectral density is shown to remove much of the broad¬ 
ening seen in the raw relaxation data in Fig.@](f), and 
thus provides a better measure of the environmental dy¬ 
namics. Some small discrepancies between the resulting 
spectrum and the theoretical result are evident, however, 


we note that whilst the theory incorporates effects such 
as hyperfine couplings and dephasing rates, more com¬ 
plicated effects such as g-factor anisotropy and enhance¬ 
ment, interactions with other paramagnetic impurities, 
and other strain-related phenomena, have not been in¬ 
cluded. Such effects are highly sample-dependent and 
thus difficult to predict in general terms; however the 
technique developed in this work provides an ideal means 
by which to facilitate their investigation. Finally, we note 
that although this demonstration involves ensembles of 
NV and PI centres, the detection is highly local as the 
response of each NV is dominated by its nearest PI cen¬ 
tre: for this sample (50 ppm PI) the mean distance to 
the nearest PI centre is about 2.7nm—. 


III. DISCUSSION 

In this work we have presented a general method for 
extracting the spectral distribution of an arbitrary elec¬ 
tronic environment based on tuning a spin -1 NV probe 
system, via controlled application of an external field, 
into resonance with the transitions of the target elec¬ 
tronic system. The method was tested using an ensem¬ 
ble of NV centres in a type-lb diamond sample to de¬ 
termine the spectral distribution of the PI spin bath, 
showing excellent agreement with the theoretical expec¬ 
tations. Our relaxation based ESR method has a number 
of advantages over existing techniques. Measurements of 
the NV relaxation in general do not require microwave 
control, and thus require no manipulation of the sam¬ 
ple. With relaxation times much longer than dephasing 
times, Tj based protocols can be significantly more sen¬ 
sitive to ESR detection. Finally, even in the ensemble 
case demonstrated here, the NV spin relaxation is domi¬ 
nated by local interactions with the environment affords 
an effective spatial resolution of a few nanometres. By 
extending this technique to the single-NV probe case (or 
other solid-state single spin systems), determination and 
characterisation of nanoscale ESR spectra of single elec¬ 
tronic systems will be possible. 
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Appendix A: Modeling of the NV relaxation process 


Here we develop the theory of relaxation based sensing as used in the main text of this work. As we are considering 
axial magnetic field strengths, Bq, such that Bq ~ D/2 ~ 512 G, only the |0) -SA | — 1} transitions of the NV spin 
will be appreciably excited, meaning we can disregard any population of the | + 1} state. The time evolution of the 
associated density matrix is described by 


^ = -i[n T (t), PT \, (Ai) 

where pt represents the combined density matrix of the entire spin + environment system. The full Hamiltonian 
is given by Ht = Bq + V + He where Hq and He are the self Hamiltonians of the NV centre and environment 
respectively. The coupling of the environment to the NV is described by the full dipolar interaction due to all spins 
in the environment: 


V = 


Po_ 

4-7T 


hjnvlE ^3 ^ _ 3 ^2 ' R *) ( R * ' «%) 


(A2) 


which includes both transverse and longitudinal components, proportional to 7* x ,y and V z of the NV spin respectively. 
The latter have a pure dephasing effect, resulting in an additional contribution to the intrinsic dephasing rate of the 
NV. As relaxation processes occurr on timescales that are much longer that the typical interaction timescales of the 
environmental constituents, the resulting dephasing will be purely exponential. These effects may thus be modeled 
using a master equation approach for the reduced density matrix for the ensemble averaged dynamics of the NV, p, 
as follows 


+ E (ap4- 

NV,i ' 

where, in the present context, Ci is the Lindbladian operator corresponding to a pure dephasing process on spin i, and 
is given by £, = y/TT^Sz^. The total dephasing rate due to both the local crystal environment and the longitudinal 
coupling to the environment is given by T 2 = (X^)^ 1 + T^ - ®. The timescale of the intrinsic dephasing process is 
described using the inhomogeneous linewidth, (TJ) -1 , since the transverse phase accumulation occurs in the absence 
of any pulsed microwave control. Subtle tuning effects that modify the sensitivity of this technique to various parts of 
the environmental spectral density may be achieved by changing the intrinsic dephasing rate via dynamic decoupling 
techniques. 

In what follows, owing to the strong intra-environment and comparatively weak NV-environment couplings, we will 
treat the coupling of the environment to the transverse components of the NV spin as a semiclassical oscillatory field 
(these simplifications will be rigorously justified in what follows), 

V = Be iut +Bh~ iut , (A4) 

where B = B x (u>)S x + B y (ui)S y ; and B x and B y are the x and y components of the magnetic field. The frequency 
spectrum is determined by analysing the interaction between environmental constituents, as described by He- To 
make the solution tractable, we change to the interaction picture. The transformed equation of motion is given by 

+T 2 (S z p I (t)S z - pi{t)), (A5) 

with the interaction Hamiltionian given by V/ = e* Wot Ve _lWo *. 

We are then interested in determining the rate at which the NV spin relaxes to its equilibrium state under the 
influence of the environment. We proceed by reducing the 3x3 system of first order linear differential equations 
described by eauation lAhl to a higher order differential equation for Pq = poo- We then wish to solve this equation, 
together with the initial conditions of py = 0 unless i = j = 0, in which case we have poo = 1, representing the initial 
polarisation of the NV spin in the |0) state. 

To gain insight into the expected analytic solution for the spin-1 NV centre, we consider the simplified case in 
which only one of the transitions of the NV centre is excited by the environment, and the other is assumed to be 
too far detuned to have any effect on the population of the spin states. This simplifies the analysis dramatically, yet 
demonstrates the main properties of relaxation based detection. Incidentally, this simplification is still applicable to 
a spin-1 system for cases of significantly strong Zeeman splittings between the | ± 1) states. This ensures that one 
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transition will be excited by the environment, whilst the other is not. This forms the basis of the technique by which 
environmental spectra may be mapped. 

The equation of motion for P 0 ( w Ej An B,t) = poo{t) is 

+ 2T+ (T 2 + S 2 + 2P 2 ) ^ + 2TP 2 P 0 - TP 2 = 0, (A6) 

at 6 at z v ' d£ 

where B = (P 2 ) 2 is the second moment of the strength of the effective magnetic field operator, B. 


1. Response to a monochromatic transverse field 

a. Resonant case 

For the case where the frequency of the environment is resonant with the transition frequency between the probe’s 
spin states (<5 = 0), the solution of Ea. IAGI is 


, , ii / r 2 t 

Po(uj E ,Bo,B,t) = — + — exp (—— 


cosh -a/To — 8B 2 ) + _ 

' 4 V 2 ) ^r 2 - 8P 2 


r 2 


sinh - A /r 2 -8P 2 


(A7) 


Typically the spin based environments in which we are interested couple weakly to the NV spin as compared with its 
intrinsic dephasing rate, implying T 2 P. In fact, even a strong coupling will also induce additional dephasing, so 
even in a worst case scenario we are guaranteed T 2 > P. In this limit, we have 


Po(B,t )\ s= o = i ^ exp (-2^) • 


Hence the resonant (and therefore maximal) longitudinal relaxation rate is given by 

d 2 

r max = 2 - = 2 T* (B 2 ) . 

r 2 


(A8) 


(A9) 


b. General case 

When a finite detuning, 6, exists, we may examine relative importance of the terms within Ea lAfil subject to 
rescaling t in terms of the decay time from the resonant solution. That is, if we consider the dimensionless variable 

T = rf“t and retain terms up to and including order O (fl ) , the solution tor an arbitrary detuning becomes 

a>(S,t)| 5=0 = 1 + lexp f -llb'-l pAyr ) - ( A1 °) 

For zero detuning, we recover the previous result ('Ea. lA9l) . For finite detuning, the relaxation rate is modified by a 
Lorentzian factor with a FWHM of T 2 . The complete decay profile is then obtained by integrating this expression 
over the spectral density of the environment, implying that the d-dependent relaxation rate acts acts to filter out the 
environmental spectrum about S ~ 0. 


2. Response to a transverse field with an arbitrary spectral density 

Even without considering the specifics of the spectral density, the response of the NV spin to an arbitrary spin 
bath can vary remarkably due the geometric proximity and arrangement of the bath relative to the NV centre. The 
definitions of the NV spin relaxation and the corresponding filter function are given by 


M(P 0 ,t) = 


G(w E ,wo,i)5'(wE,wo)da;E, 


(All) 











10 


,and 

G(wE,wo,i) = J Po(we, wo, B, t)PB(B) dB, (A12) 

respectively, where the filter function, G(we, wo) acts t° filter out regions of the spectral density (as dependent on the 
external field strength, wo), and depends explicitly on the geometric arrangement of the environmental constituents. 
Ultimately, given some measurement record and filter function, it is expression IA11I that must be deconvolved to 
reproduce the spectral density, S(we). In this section, we consider the effects of the geometric arrangement of the 
environment on the filter function, G, for a general spectral density; and specific cases of important spectral densities, 
namely that due to the internal PI nitrogen donor electron spin bath in type-lb diamond and its corresponding 
surface spin bath, are considered below in SectionlCl 


a. Response to an internal (bulk) spin bath 


We note that the coupling of the NV to a bath spin located at some distance r may be written B = b(9, (f>)/r 3 , 
where the specific details of the angular dependence of the coupling are incorporated into the parameter b(9, <f>). We 
note that is is necessary to omit further discussion of b until SectionlCl as different environmental processes will be 
more readily detectable by the NV spin at different relative angles. Unlike the transverse (spin echo) case, where 
the precession of the NV spin vector in the x — y plane is sensitive to all longitudinal held sources, the effect on 
the longitudinal projection is dominated by the coupling to the nearest PI electron spin. As such, we may write 
Pb{B) = P r (r)Pe.d>(0, <t>), where are the spherical coordinates associated with the distribution of held sources. 

From Ref.l29l the distribution of distances from a given NV centre to its nearest spin impurity is given by 


where n is the 
M(u 0 ,t) = 


P r (r) = 47rnr 2 exp ^— 7 rnr 3 ^ , 


(A13) 


average density of impurities in the bath. Substituting this expression into Ea. IAllI we hnd 
1 f 


-V2irb(9, <j))n 


I r 2* ^ 3,0 f 8& 2 (6>,</>)n 2 7r 2 fr 2 

n + 52 ^0,3 


h2 (9, (j>)n 2 TT 2 t 

9 (r| + <5 2 ) 


-±,0,i \ P e ^(9,(j))d9d(/)S(u E ,ujo)dLo E 


/ exp (~T ( 27rn2f ) 1/2 (l b D y r^ s 2 ) 5 ( WE ’ w °) d We , 


(A14) 


for smal t, where 0 is the Meijer-G function, and (|6|) = / \b(9, <f>)\Pg^(9, <j>) d9 d<j>. Thus, we identify the hlter 
function associated with environments inside the diamond lattice to be 


Gin(<jJE,Wo) = 


(we — d wo) 2 + r! 


(A15) 


where Ai n is a constant, associated with the geometry of the bath, that may be renormalised. 


b. Response to an external (surface) spin bath 

In contrast to the bulk spin bath case, spins on the surface are unable to exist arbitrarily close to the NV centre. 
Typically we consider samples in which NV centres exist at some depth h + Sh below the surface, with h being the 
mean depth, and Sh a normally distributed variable with variance (Sh 2S ) C ft 2 . In this case, an individual NV spin is 
exposed to many bath spins, meaning that the effective coupling distribution, Pb(-B), is normally distributed. 

In this case, we may expand Eci. lAlOl for small t, which, upon substitution into Ea. IAllI and averaging over P b(B) 
gives 

(Po(t)) ~ l-dt(B 2 )(j f ^^5(w E ,wo)dw E ) (A16) 

Thus, we identify the filter function associated with environments outside the diamond lattice to be 

r 2 _ 

(we — D + co 0 ) + 


G 0 ut(wE,Wg) — A ont 


(A17) 
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Appendix B: Derivation of the deconvolution algorithm 


Although we wish to invert Ea. IAlll for the spectrum So = S{lu e, 0), we must also deal with the intrinsic shot noise, 
r)(Qo) that arises during the measurement process. Thus, we are interested in solving 

M(fi 0 ) = {So * G) (f2 0 ) + h(^o)) (Bl) 

where the effective spectral features of may be inferred from the fact that the errors associated with any two 
measured time points are uncorrelated from one-another. Our goal is thus to find some function, H , such that a 
least-squares estimate of So is given by 

S 0 = H * M. (B2) 


In order to deal with products of the relevant functions, rather than convolutions, we switch to Fourier space via the 
Fourier transform, denoted 1?. 

The mean-square error in So is given by 

LSE = Ex{|3(5 0 ) -S(So)| 2 } 

= Ex{|3fa03(ri)-3(s 0 )| 2 } 

= Ex {IS (H) [3 (So) 3 (G) + 3 ( r ,)] - 5 (S 0 ) | 2 } 

= Ex { (3 (So) [3 (H) 3 (G) - 1]) (3 (So) [3 (H) 3 (G) - 1])*} + Ex { (3 (H) 3 fa)) (3 {H) 3 fa)) *} 

+Ex {(3 (So) [3 {H) 3 (G) - 1]) (3 {H) 3 fa)) * } + Ex {(3 {H) 3 fa)) (3 (S 0 ) [3 {H) 3 (G) - 1])* } 

= [3 {H) 3 (G) - 1] [3 {H) 3 (G) - 1]* Ex { 13 (So) I 2 }+ 3 {H) 3* (H) Ex { |3 fa) | 2 } , (B3) 

where the last line follows from the fact that the noise and the environmental spectrum are uncorrelated. Minimising 
the least-squares estimate of this quantity with respect to 3 (H) (see Section[B]), we find 


3*(G) 13(S 0 )1 2 

|3(G)| 2 |3(S 0 )| 2 + |3fa)| 2 ' 


(B4) 


Thus, we have 


So = 3 _1 (3(H)3(M)). 


(B5) 


This shows that determination of the appropriate deconvolution filter requires knowledge of So, which is what we are 
actually trying to measure, however, since the measurement record, M{ui o) contains many spectral features similar to 
So, M{u o) provides an ideal starting point with which to implement an iterative procedure using this method. 


Appendix C: Theoretical description of the coupled NV-P1 system 

In this section, we discuss the features we expect to be evident in the PI centre spectrum by examining the effect 
of a PI centre on the magnetic field dependent relaxation rate of a near-by NV centre. We conclude this section by 
demonstrating the equivalence of the semi-classical approach used in this work, and a quantum mechanical treatment 
of the NV-P1 interaction. 


1. PI Hamiltonian 


The Hamiltonian of a PI centre is given by 

n P1 = s-A P1 -i +B 0 - ( 7 eS + 7 N l) +f-Qpi-f, (Cl) 

where Api is the hyperfine tensor describing the coupling between the PI electron, S , and 14 N nuclear spin, I; and Q 
is the quadrupole splitting of the nuclear spin. For all field strengths at which there is an appreciable overlap of this 
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spectrum with the NV spin filter function, Bq ~ 512 G, we find that the eigenstates of the PI centre electron spin are 
predominantly dictated by the external magnetic field. In this instance, the Hamiltonian of the PI centre becomes 

H$i = A z S z l z + A x {S x l x +S y l y ) + B 0 ( le S z + lN l z ) + Q P1 l 2 z , (C2) 

for cases where the PI axis is aligned with the external magnetic field, where A z = 114MHz, and A x = 81.3 MHz. If 
the PI axis is aligned along one of the three other bond directions not aligned with the field, the Hamiltonian may 
be transformed via the rotation operator 1Z = exp (—il v 9) = 1 — il y sin(0) — I' 2 (1 — cos(0)) (The other two axes may 
be realised via a trivial rotation about the z axis), and is given by 

T~Lp i = O'z^z^-z + a x {S X I X + Syly) + .Bo (%<S Z + 7n^z) + (C3) 

where a z = | (8H X . + Az) = 85 MHz, and a x = | (5 A x + 4 Az) = 99 MHz. 


2. Coupling of the PI environment to the NV centre 


The interaction between the NV spin and the PI nuclear spin is ignored on account of its comparative weakness. 
For the interaction between the NV and the PI electron, the magnetic dipole interaction is given by 


'Hint = B u 


P'5-3 r-P f S 


where r is the unit separation vector between the NV and PI centres, and 


^int 


47rr 3 ’ 


(C4) 


(C5) 


is the effective dipolar coupling strength. 

The components of the NV-P1 interaction responsible for the relaxation of the NV spin are those coupling to its 
transverse components, namely S x and S y . Without loss of generality, we may rewrite this interaction as an effective 
quantum mechanical magnetic field, B , where B XiVtZ couples to S x , y , z - 

To make use of Ea lAlll we make the semi-classical approximation, and assume that NV-P1 interaction plays very 
little part in determining the dynamics of the environment. The problem of determining the environmental spectrum 
then reduces to solving for the environmental evolution exclusively under its own influence as follows. 

The effective field strength associated with the allowed NV-P1 |0,4_) O | + 1, f) transition is given by 


-Ball = ^ \/ ( B XX + Byy) + (2 B X y) 

= ^B sin 2 ( 0 ), (C 6 ) 

and that due to the disallowed 10 , ms, mi) -O- | + 1 , ms + 1 , mi — 1 ), | 0 , ms, mi) -O- | + 1 , ms — 1 , mi + 1 ) transitions 
is given by 

Bdis = 2 \JB xz + By Z 

= ^B sin (20) (C7) 

To determine the effect of the PI environment, we compute the autocorrelation functions associated with the field 
components above. Interactions between environmental components may be modeled by damping these autocorrela¬ 
tion functions with a decaying exponential, exp (—Tpit) to describe their relaxation due to mutual flip-flop processes 
with corresponding relaxation rate Tpi. 

The corresponding spectra may then be found by computing the Fourier transforms of the autocorrelation functions. 
From this, we find the spectra associated with the allowed transitions to be 


and 


s a °5M = 


67T 


PI 


PI 


+ 


PI 


Tpr + (w ± w 0 ) 2 rp! + (w ± UJ 0 + A z ) 2 ' r| x + {to ± too - A z ) 2 


(C8) 


S a °nM = 


Tpi 


Ppi 


Tpi 


671- [Pp! + (w ± w 0 ) 2 Tpj + (oj ± oj 0 + a z ) 2 r 2 ! + (w ± w 0 — a z ) 2 _ 


(C9) 
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for cases of on and off axis PI centres respectively. Taking the relative proportions of on and off-axis PI centres to 
be 25% and 75% respectively, we find the overall spectrum associated with the allowed transitions to be 


S'aii(w) 



(CIO) 


Similarly, the spectra associated with the disallowed transitions are given by 


«<") = i 


Tpi 


Tpi 


+ (bJ zh CJq H“ Ai ) 2 ^pi + (cJ =t CJq — A2) 2 


and 


= 


47T 


Tpi 


Tpi 


Tp-j^ + (w A UJq + A 3) 2 + (w i UJg — A 4) 2 


respectively, where 


(Cll) 


(C12) 



The spectrum associated with the disallowed transitions is then 


(C13) 


By employing the full spectrum, S(oje ) = SWi^-e) + S^s(cje), in eauation lAlll we find the resulting external 
field-dependent relaxation rate of the NV centre to be 


ri(w 0 ) — 


(Bj) 

6ir 


r 2 + Tpi 


+ 


r 2 + Tpi 


Bi 


4 V ( r 2 + Tpi) 2 + 4(w 0 ± f + 4-) 2 (r 2 + Tpi) 2 + 4(wo ± f - 4") 2 

3 ( r 2 + rpi , r 2 + rpi 


+ 


4 V( r 2 + r Pi) 2 + 4 (^0 ± f + ) 2 (r 2 + Tpi ) 2 + 4(wo ± y ~ if) 2 , 

r 2 + Ppi 


(r 2 + rpi) +4(w 0 ±§) 2 


47T 


r 2 + Tpi 


r 2 + Tpi 


4 V( r 2 + r P i) 2 + 4( Wo ±f + 4^) 2 ' (r 2 + r P i) 2 + 4(u;o±f - , 

3 ( r 2 + Tpi r 2 + Tpi \ 


+ 


4 V ( r 2 + r P i) 2 + 4( Wo ± § + ^) 2 (r 2 + r P1 ) 2 + 4(cu 0 ± § - ^) 2 


(C14) 


where the effective couplings, (.B 2 ) and are due to integration over all possible NV-P1 separations. By taking 

the average PI density to be 50ppm, and the FID rate to be T 2 = 5.0MHz, we may plot the resulting field-dependent 
relaxation rate of the NV spin (Fig.[SJ). _ 
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Background magnetic field strength [G] 

FIG. 5: Analytic calculation of the relaxation rate of an NV centre spin placed in a 50ppm environment of PI donor electron 
spins vs the strength of an external magnetic field aligned along the NV axis, (a) Peaks associated with allowed transitions 
of the NV centre spin due to its direct interaction with aligned PI centre electron spins, (b) As in (a), but for the case of PI 
centres not aligned along the NV axis, (c) Peaks associated with disallowed transitions of the NV spin as mediated by flip-flop 
dynamics between the on-axis PI electron and nuclear spins, (d) As in (c), but for the case of off-axis PI centres. 










